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-^A  magnetic  electron  spectrometer  for  electrons  to  10  MeV  energy 
has  been  designed,  constructed,  and  tested.  The  magnet  assembly 
uses  ra re -ea rth/ cobalt  magnets  and  achieves  a  gap  field  of  4  -  4.  S  kG. 
The  operation  of  the  spectrometer  was  tested  with  beta  sources  and  po¬ 
sition  sensitive  detectors.  It  was  found  that  a  coincidence  detector  in 
the  entrance  collimator  eliminates  most  of  the  bremi strahlung  backgro 
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\  iO.  Abstract  (continued) 

^4ome  internal  baffles  are  also  necessary  to  reduce  scattered  electron 
background.  The  tests  showed  that  a  magnetic  spectrometer  which  can 
analyze  electrons  to  more  than  10  MeV  can  be  readily  constructed.  A 
preliminary  design  for  a  flight  unit  magnetic  electron  spectrometer  is 
given. 
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1.  INTRODUCTION 


The  reliable  detection  of  high  energy  (1  to  10  MeV)  electrons  in  space 
has  been  a  difficult  undertaking.  Electrons  scatter  considerably  in  passing 
through  matter  and  produce  bremsstrahlung  which  can  interact  to  simulate 
electrons  of  different  energy.  Most  regions  of  the  earth' s  magnetosphere 
which  can  contain  high  energy  electrons  also  contain  high  energy  protons, 
and  the  electrons  become  difficult  to  detect  in  the  presence  of  the  other  rad¬ 
iations. 

A  number  of  approaches  have  been  used  in  the  past  to  detect  high  ener¬ 
gy  electrons  from  satellites.  One  of  the  first  such  detectors  was  a  magnet- 
scintillator  combination  for  0.05  to  4  MeV  electrons,  which  was  used  along 
with  a  set  of  absorbe  r-dete  rnuned  threshold  detectors  for  1  to  5  MeV  elec¬ 
trons  (Re:.  1.  1).  The  magnet/  scintillator  combination  was  used  to  be  able 
to  iistinguish  electrons  from  protons,  but  bremsstrahlung  induced  back¬ 
ground  produced  most  of  the  counts  below  1.5  MeV  (Ref.  1.2).  The  absor¬ 
ber-determined  threshold  detectors  were  also  sensitive  to  protons  above' 

20  to  35  MeV  and  so  cannot  distinguish  1  to  5  MeV  electrons  from  these 
higher  energy  protons. 

A  l  to  10  MeV  electron  spectrometer  has  been  constructed  using  a 
dE/ dx  -  scintillator  coincidence  measurement  (Ref.  1.3).  This  combin¬ 
ation  is  sensitive  to  high  energy  protons  (  >  70  MeV)  and  so  has  limited  use 
in  much  ot  the  magnetosphere.  A  magnetic  spectrometer  for  0.  1  to  1  MeV 
electrons  has  also  been  constructed  (Rcl.  1.4).  The  magnetic  spectrometer 
used  a  gap  field  of  about  1  k  Gauss  and  so  was  limited  to  about  1  MeV  elec¬ 
tron  energy  to  avoid  excessive  size  and  weight.  The  magnetic  spectrometer 
had  no  c  oincidence  detector  and,  therefore,  is  susceptible  to  penetrating 
high  energy  protons  and  bremsstrahlung. 

A  magnetic  electron  spectrometer  for  energies  up  to  2.5  MeV  has  been 
flown  (Ref.  1.5)  as  has  one  for  energies  up  to  3.  1  MeV  (Ref.  1.6).  Both  of 
these  magnetic  analyzers  use  180  degree  bending  and  have  no  coincidence  de¬ 
tector.  They  are  thus  susceptible  to  high  energy  penetrating  protons  and 
to  bremsstrahlung  background.  The  magnetic  spectrometer  in  Ref.  1.5  is 
reported  to  have  low  internal  bremsstrahlung  sensitivity  because  of  exten¬ 
sive  baffles,  but  external  bremsstrahlung  flux  levels  can  be  quite  intense, 
require  large  amounts  of  lead  (high  Z)  shielding,  and  become  much  more  se¬ 
vere  at  higher  electron  energies.  The  magnetic  spectrometer  in  Ref.  1.6 
used  shielded  detectors  adjacent  to  unshielded  detectors  to  measure  back¬ 
ground.  This  leads  to  gaps  in  the  measured  spectrum,  as  well  as  limiting  the 
measureable  electron  count  rates  to  a  fraction  of  the  background  count  rates; 
the  actual  fraction  depending  on  counting  statistics  and  background  uniform¬ 
ity. 

Advances  in  magnetic  materials  make  it  now  possible  to  construct  mag¬ 
netic  analyzers  with  much  larger  gap  fields,  up  to  nearly  5  kG.  It  thus  be¬ 
comes  possible  to  construct  an  instrument  to  detect  up  to  10  MeV  electrons, 
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and  which  is  not  much  larger  than  previous  designs.  There  are  possible 
simplilications  from  using  position  sensitive  detectors;  therefore,  a 
thorough  investigation  of  the  design  of  a  1  to  10  MeV  magnetic  electron 
spectrometer  is  desircable. 

This  report  describes  the  design,  construction,  and  testing  ol  sev¬ 
eral  <  ratio!.:.,  i,.r  .»  magnetic  electron  ip<  ter  to  cover  to  at 

least  10  MeV.  The  results  re  used  to  develop  a  preliminary  design  of  a 
flight  unit. 
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2.  T UK  MAGNET  ASSEMBLY 


V 


Z.  1  Design  and  Construction 

A  general  dc -.lgn  for  an  electron  magnetic  analyzer  is  shown  in  Figure 
Z.  1.  This  design  includes  a  coincidence  detector,  a  sweeping  magnet  for  low 
energy  electrons,  and  a  proton  telescope  located  in  a  position  completely  shielded 
i  rom  electrons  to  ai  least  1'  MeV.  I  he  design  uses  a  semi-circular  array  of 
n  ignets  oil  two  iron  plates,  with  a  completely  enclosing  outside  return  yoke. 

The  actual  magnet  design  for  testing  purposes  is  shown  in  Figure  2.2. 

I  his  show  s  the  outline  ol  the  top  and  bottom  iron  plates,  the  edge  pieces  which 
provide  support  and  a  magnetic  return  path,  and  the  actual  positions  of  the  Sm- 
Co  magnets.  The  iron  plates  and  iron  edge  pieces  are  made  of  1/4  in.  h  thick 

,  which  it  a  hi  penm  ty,  ft  iron.  Becaux  •  <  |  ** 

iron,  the  attractive  force  ol  the  two  ■  agnet  poles,  a  set  of  aluminum  support  bars 
is  bolted  to  each  plate  set  (composed  of  two  three-inch  wide  pieces  :  l/4  inch 
•.n,n.,  a  «  .  .  '■  '•••  .*\  -■  I  <  w 

for  *tt.u  hment  to  a  fixture  used  to  ring  the  two  magnet  poles  together  in  the  fi¬ 
nal  magnet  assembly. 

The  ; .  ignrt  design  uses  lh  magnets  of  nominally  1  x  1  \  1/2  niche  and 
2  ]  Inchi  ■  <  ■  ■  A 

tb  sq  i  . .  •  r:>  re-earth  cobalt  ceramic  mac  nets  (  1  l  1 1  x  1  1/  '  b  x  l/2  in*  hes 
wa  n >  d  !  r  the  p*  los.  The  magnets  win  first  measured  for  field  strength  at 
j,v,  ..  si ti  '  o*  e  .  h  side  of  the  magnet,  and  sorted  int  *  two  groups,  based 
:i  t'  e  ;hes|  vor.ige  field  't  rengt  .  f*  ’  ■  »  >p  and  bottom  pole  i.  These  mag* 

•  .  ;  r *  t  nrn  a  r ranged  t*>  give  th  *  ■  t  r *  r  _  e  s  t  .  i.*  ar  thr  outer  edgi*  1  the  ma i_  - 

r  t  pole,  where  the  highest  energ-.  (10  MeV)  elr  tr  u.s  travel,  to  ensure  that  the 
magnet  assembly  would  analyse  t  •  at  least  10  MeV . 

Placing  the  magnets  onto  th<  pole  plates  requires  considt  rable  *  are.  I  he 
magnets  are  strongly  altra*  ted  to  thi  ir  >n,  .«nd  >erause  of  their  parallel  magneti- 

.  ,  ■  |  ,  1 1  I  I  .  ■  I  1  Bet  ■  :  1  -•!••*  w  '  .i!  brittle, 

that  a  s<  ha  2  eel  •  tipping*  Each  pole  is  as  tern - 

Ted,  nr  magnet  at  a  time,  slipped  or.to  the  iron  plate  near  the  edge  to  avoid 
itrong  shocks.  The  outer  edges  of  the*  magnets  are  held  in  place  by  aluminum 
nracket  and  stainless  steel  screw  A**  the  magnets  are  added  toon*  el  c,  screw 
and  washers  are  used  to  hold  them  in  place  against  the  mutual  repulsive  forces. 

;  ha  arrangement  ol  the  internal  screws  is  shown  in  Figure  2.2,  wuh  all  oi  the 
outer  edges  of  the  magnets  being  ;.<Td  by  aluminum  brackets.  The*  magnet  arrange- 
m,  2. 2  has  spaces  between  the  threi  rows  of  magnate  because  of  the 

necessity  of  using  screws  to  hold  the  magnets  against  their  mutual  repulsion 
during  assembly  of  the  individual  poles. 

On.  ••  the  two  poles  are  completed  with  the  magnets,  the  1  1/2  inch  high  iron 
edge  piec  es  are  mounted  on  one  pole.  The  poles  are  then  bolted,  via  the  alumi¬ 
num  support  pieces,  to  the  plates  of  an  assembly  rig,  as  illustrated  in  Figure  2.4. 
The  pole  pieces  are  bolted  in  place  with  the  assembly  rig  taken  apart,  and  the 
top  plate/ moveable  plate/ slide  rods  component  is  then  placed  on  top  of  the  bot¬ 
tom  plate.  The  magnet  poles  are  separated  by  at  least  five  inches  during  this  step, 
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Figure  2.  1  General  Design  of  a  Magnet  Assembly  for  the  Electron 
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Fin.  1.  4.  Assembly  Rig  for  Magnet  Pole  Pieces 


since  otherwise  the  magnet  poles  attract  with  too  much  force  to  allow  the 
operation  to  be  easy  and  safe.  After  the  bottom  plate  has  been  securely 
fastened,  the  magnet  poles  are  slowly  brought  together  by  turning  the  nut 
on  the  threaded  rod.  The  assembly  rig  must  be  rigid  and  strong,  since  the 
magnet  poles  attract  with  nearby  1000  lbs.  of  force  just  before  contact  of 
the  top  pole  piece  and  the  edge  pieces.  After  the  top  pole  has  been  fastened 
to  the  edge  pieces  with  screws,  the  complete  magnet  assembly  is  removed 
from  the  assembly  rig  and  is  ready  for  testing. 

2.2  Magnetic  Field  Properties 

The  completed  magnet  assembly  has  a  gap  of  1  /  2  inch  between  the  mag¬ 
net  faces  and  a  clear  gap  of  about  3/8  inch  between  the  aluminum  magnet  edge 
brackets.  With  an  iron  plate  in  the  detector  insertion  opening  in  the  top  plate 
(Figure  2.  3),  the  two  removeable  edge  pieces  (Figure  2.2)  were  taken  out  one 
at  a  time  as  necessary,  and  the  magnetic  field  profile  measured  over  the  entire 
magnet  pole  area  in  the  center  of  the  magnet  gap.  The  resulting  field  profile 
is  shown  in  Figure  2.5. 

The  magnetic  field  profile  has  a  number  of  characteristics.  The  field 
in  the  center  of  the  outer  semicircle  of  magnets  averages  slightly  over  4.5  kg 
and  shows  that  this  field  intensity  can  be  readily  achieved  with  the  present  de¬ 
sign.  The  lower  field  intensity  at  the  inner  part  of  the  semicircle  is  because 
the  weaker  magnets  were  placed  there.  By  selection  from  a  larger  number 
of  magnets,  it  should  be  possible  to  bring  the  field  at  the  center  of  all  magnets 
to  at  least  4.5  kg.  This  was  not  done  with  the  test  magnet  assembly  because 
no  extra  lxl  inch  magnets  were  available. 

The  junction  line  between  adjacent,  touching  magnets  shows  a  drop  in 
field  intensity  of  5  -  10%.  This  shows  up  at  nearly  all  magnet  junctions  and  in¬ 
dicates  that  even  with  magnets  selected  from  a  larger  number,  there  will  still 
be  variations  of  up  to  5%  from  the  average  field.  Of  greater  importance  are 
the  more  than  15%  variations  that  show  up  where  the  magnet  holding  screws 
have  forced  a  gap  between  magnet  rows.  It  is  expected  that  this  can  be  reduced 
by  having  the  corners  on  the  magnets  ground  slightly  so  that  the  holding  screws 
can  be  placed  at  the  corners  where  four  (4)  magnets  meet,  and  so  eliminate  all 
gaps  at  magnet  sides.  By  this  design  change,  a  flight  magnet  assembly  can 
be  produced  w-ith  5%  field  uniformity  over  the  entire  magnet  semicircle. 

The  field  variations  of  5%  are  not  expected  to  affect  significantly  the  per¬ 
formances  of  the  magnet  assembly  for  electron  energy  analysis.  For  a  5  MrV 
electron,  field  variations  of  5%  over  a  1  cm  distance  produce  only  a  0  1%  change 
in  deflection  radius,  and  this  is  safely  neglected.  The  change  in  exit  distance 
for  1  MeV  electron  is  still  only  about  1  -  2%,  which  is  smaller  than  the  fringing 
field  corrections.  Since  it  is  expected  that  a  flight  unit  would  have  to  be  cali¬ 
brated  with  electron  sources  and  beams,  the  field  variations  of  5%  should 
have  no  net  effect  on  the  performance  of  a  flight  instrument. 

The  field  profile  along  the  particle  entrance  line  (the  entering  electron  line 
of  Fig.  2.2  extrapolated  straight  through  the  magnet  gap)  is  shown  in  Figure  2.6. 
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Measured  Magnetic  Field  Intensities  of  the  Assembled  Analyzing  Magnet 


This  prolilc  shows  the  20%  field  dip  in  the  1/  16  inch  gap  between  the  magnet 
rows,  which  should  be  reduced  to  about  6%  by  the  magnet  modification  dis¬ 
cussed  above.  Because  of  the  large  (4.  5  cm)  overhang  of  the  iron  pole  plates, 
the  fringing  field  cuts  off  more  sharply  and  actually  reverses  direction  at 
about  +1.4  cm  (0.0  cm  coincides  approximately  with  the  outer  edge  of  the  first 
magnet).  This  sharp  cut-off  and  reversal  somewhat  reduces  the  fringing  field 
effects  on  incoming  particles. 

The  measured  profile  in  figure  2.  b  ha*  been  fitted  to  the  first  maximum 
at  -2.0  cm  and  used  in  fringing  field  calculations.  The  procedure  is  discusse 
in  detail  in  Appendix  A.  For  precise  calculations  of  electron  paths  in  the  mag¬ 
net  assembly,  both  the  f ringing  tie'!  effects  and  the  variation  of  the  magnetic 
field  over  the  pole  face  (Figure  2.  5)  must  be  taken  into  account.  This  is  par¬ 
ticularly  true  for  a  flight  magnet  .  .sembly,  where  onl\  5°c  field  variations  are 
expected,  for  the  test  magnet  a*  *mbly,  wher-  the  Held  variations  are  20’', 
and  a  lower  intensity  field  occur  near  much  of  the  low  energy  cle  tron  exit 
portion,  the  major  effect  on  electron  paths  conies  trom  the  interna  field  var¬ 
iations. 

The  approximate  paths  an  exit  positions  for  1  -  10  MeV  <  c-  ti  :  ..  r> 

shown  in  .  tgure  2.  7.  The  elec  tron  paths  were  calc  elated  using  ivi  :.-si 

field  intensity  along  the  semicircular  part  of  cad  ath  be-  en  ia.m~.et  poi  s. 
The  mere  .se  m  Held  intensity  on  the  cuter  semicircular  p  t  oi  the  :nu.  net 
pole  results  i<  the  spacing  change  it  fc  MeV.  A*  shown,  <'lcctr  ns  of  exactly 
10  MeV  energy  will  likely  not  b>  analysed  but  will  c  scap  near  the  back  edge 
thc  pole  piece.  This  re  sults  primarily  from  the  reduction  in  field  intensity  a: 
the  1/ lb  inch  gaps  between  mag:.et  r<  a  s  (Figure  2.2)  and  snou.d  not  ccur  in 
a  flight  magnet  assembly  w  hi  *  h  J  •  •  not  ha-.  <  th<  s<-  gaps.  1  he  i  mil;  n  cl  a 
tew  extra  lxl/2  and  1/2  x  1/2  inch  magnets  would  also  make  In  magret  as¬ 
sembly  capable  of  analyzing  to  at  le  ast  10  MeV  electrons.  The  test  magnet 
assembly  was  only  used  with  beta  (electr  .n)sources  to  3.5  MeV,  so  the  ques¬ 
tionable  analysis  at  10  MeV  is  no  problem.  Since  the  test  magnet  assembly 
definitely  analyzes  to  >9  MeV,  the  effect  of  lower  energy  electron  and  brrm- 
sstrahlung  background  on  the  analysis  ol  electrons  near  10  MeV  can  still  he  in¬ 
vestigated. 

Figure  2.7  also  shows  the  outline  of  part  of  the-  entrance  collima  and 
the  location  of  the  position  sensitive  detector  in  the  lowest  energy  position. 

For  this  test  set-up,  the  position  sensitive  detector  has  a  lowest  detectable 
energy  of  about  1.6  MeV,  and  a  maximum  of  about  4.4  MeV.  All  ol  the  test¬ 
ing  was  done  with  a  Sr-  ’i  -  90  beta  source  of  2.27  MeV  maximum  electron  energy 
and  a  Ru-Rh-lOb  beta  source  of  3.54  MeV  maximum  electron  energy.  Addi¬ 
tional  details  of  the  test  set-up  and  the  results  of  the  testing  are  given  in  the 
following  section. 

3.  THE  MAGNETIC  ELECTRON  SPECTROMETER 
3.  1  Spectrometer  Assembly  for  Laboratory  Tests 

The  tost  magnet  assembly  described  in  Section  2  was  used  with  position 
sensitive  detectors,  a  200  pan  totally  depleted  surface  harrier  detector  for  coin- 
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Figure  2.7  Approximate  Paths  and  Exit  Positions  for  1  -  10  MeV  Electrons  in  the  Test 

Magnet  Assembly 


cido  nee,  a  collimator  assembly,  and  internal  anti  -  scatte  ring  baffles,  as 
shown  in  Figure  3.  1.  The  test  set-up  in  Figure  3.  1  shows  the  final  config¬ 
uration  arrived  .»t  through  many  experimental  tests  with  intermediate  config¬ 
urations.  The  position  sensitive  detector  (PSD  1)  (8  nun  x  47  mm  sensitive 
area)  was  used  in  any  one  of  the  three  (3)  positions  shown,  giving  a  slight 
overlap  between  adjacent  positions.  A  second  (8  mm  x  27  mm)  position  sen¬ 
sitive  detector  (PSD  2)  was  also  used,  both  alone  and  in  anticoincidence  with 
PSD  1.  When  both  PSD’  S  were  used,  PSD  2  was  used  only  as  a  detector  since 
only  one  position  sensitive  detector  analyzer  (PSDA)  module  was  available. 

A  block  diagram  of  the  electronics  used  with  the  test  set-up  of  Figure 
3.  1  is  shown  in  Figure  3.2.  The  complete  set-up  was  not  always  used,  with 
the  200  pin  detector  part  of  the  electronics  chain  and  the  MCA  coincidence 
input  not  being  used  for  some  ol  the  earlier  tests.  When  PSD  2  was  used  in 
anti -coincidence  with  PSP  1,  it  was  connected  in  place  of  the  200  pm  detector 
and  the  MCA  run  on  anti -coincidence. 

3.2  Spectrometer  Tests  with  Radioactive  Sources 

The  experimental  configuration  for  most  of  the  beta  source  tests  with 
the  spectrometer  is  shown  in  3  igure  3.  3.  One  of  the  first  tests  was  done  with 
the  d  x  27  mm  PSD  (PSD  2)  used  with  the  PSDA,  and  positioned  to  detect  1.  5  to 
3.2  MeV  electrons.  The  resulting  spectra  (energy  and  position)  for  a  Sr-Y-90 
beta  source  are  shown  in  Figure  3.4.  The  PSD  spectrum  covers  channels  1  t«^ 

80  (128  is  lull  scale),  so  it  is  clear  that  electrons  art'  being  detected  well  be¬ 
yond  the  2.27  MeV  end-point  of  the  Sr-Y-90  beta  source.  The  log-counts  scale 
shows  that  changing  the  energy  threshold  for  position  calculation  from  250  to 
500  keV  lowers  the  electron  count  rate  at  1.5-2  MeV  log  a  factor  of  5,  while 
at  3  MeV,  it  is  lowered  only  a  factor  of  2  -  3.  The  high  energy  counts  are  a 
combination  of  bremsstrahlung  background  and  scattered  electrons  and  have  a 
higher  average  energy  loss  than  real  electrons  in  the  2  MeV'  region.  These 
tests  were  made  without  the  200  pm  detector  or  the  baffles. 

The  above  tests  were  repeated  with  the  8  x  47  mm  PSD  (PSD  1)  set  to 
detect  1.7  to  4.  5  MeV,  and  position  spectra  taken  with  Sr-Y-VO  (2.27  MeV  max¬ 
imum  electron  energy)  and  Ru-Rh-106  (3.54  MeV  maximum)  source.  The  re¬ 
sults  are  shown  in  F  igure  3.5,  and  simply  extend  the  results  shown  in  Figure  3.  4. 
All  spectra  show  counts  at  energies  above  the  beta  source  end-point.  The  Sr-Y-90 
spectrum,  which  peaks  at  about  0.8  MeV,  shows  a  definite  change  in  Blope  at  the 
end  point,  while  the  Ru-Rh-106  spectrum,  which  peaks  at  about  1.25  MeV,  shows 
only  a  slight  slope  change  at  the  end-point. 

Several  tests  were  made  to  determine  the  source  of  the  higher  energy 
pseudo-electron  background.  The  8  x  27  mm  PSD  2  was  mounted  in  a  higher 
energy  position,  and  spectra  of  PSD  1  was  taken  both  free  and  in  anti-coincidence 
with  PSD  2.  No  significant  changes  were  observed  indicating  that  electron  back- 
scattering  from  the  PSD  to  a  higher  energy  position  (>2  MeV  higher)  was  not 
the  major  source  of  background.  Increasing  the  collimator  opening  diameter  in 
front  of  the  source  by  a  factor  of  2  changed  the  position  spectrum  intensity  by 
a  factor  of  about  4,  but  did  not  change  the  shape.  This  indicates  that  much  of 
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Fig.  3.  3.  Experimental  Configuration  for  Beta  Source  Tests. 


rKry  Spectra  (left )  and  Posit  ion  Spectra  (r  ight )  for  a  lOmCiSr-Y 
ce.  for  the  appr  oxirr.at  e  posit  ion-  enei «y  range  of  1.5  to  3.2  MeV 


the  hi  h  •  mT(!y  background  is  associated  with  the  entering  electron  flux  and 
less  with  the  external  bremsstrahlung  background. 

A  set  ol  six  (6)  anti- scattering  strips  were  placed  against  each  of  the 
magnet  poles.  The  strips  were  1/  16  inch  thick  and  1/8  inch  wide,  made  of 
aluminum,  and  placed  in  a  lan  arrangement  radiating  from  slightly  behind 
PSD  1.  These  gave  no  measureable  change  in  the  position  spectra,  and  so  were 
not  effective  in  reducing  electron  scattering  and  bremsstrahlung  production. 

A  complete  set  of  position  sp«ctra  for  the  Sr  Y-90  and  Ru-Rh-106  souri  •  s 
is  shown  in  Figure  3.6.  The  background  appears  to  be  composed  of  two  com¬ 
ponent*..  A  uniform  background,  evident  above  4-5  NlrV,  is  most  likely  du« 
to  internal  and  external  bremsstrahlung  and  is  very  difficult  to  shield  without 
using  large  masses  of  lead.  The  other  background  component  appears  to  be 
from  electrons  scattered  in  the  magnet  assembly  and  is  most  noticeable  in 
the  1  MeV  region  above  the  beta  spectrum  end-point.  T*h»-  spcct.-.  in  figure 
3,0  were  taken  with  a  threshold  in  the  energy  pulse  of  125  keV.  The  energy 
spectra  for  the  1.6  -  4.4  MeV  position  of  PSD  1  are  shown  in  Figure  3.7,  with 
the  Sr-Y-90  results  being  shown  for  both  free  and  in  anti-coincic  r  i  with 
PSD  2  in  a  higher  energy  position,  and  the  R  i-Rh-106  results  being  shown  only 
for  the  anti-coincidence.  The  lack  of  .  hange  in  f  re«*/ anti -coinc  id<  nee  spec 
tra  is  in  agreement  with  earlier  results  and  conclusions  that  electrons  hack 
scattered  from  the  PSD  to  a  higher  energy  p>  sition  ("»2  Mr\  higher)  are  not 
an  important  source  of  background.  It  should  be  noted  that  the  MCA  channel 
for  0”  volts  input  was  different  n  r  the  data  in  figures  3.6  and  3.7  ("0"  V 
Channel  No.  10)  than  for  the  earlier  data  in  t  igures  3.4  and  3.  ("0"  V  =  Chan¬ 

nel  No.  0),  and  this  accounts  lor  the  peaking  at  the  lowest  position  signal  in 
k igure  3.  6. 

Ihc  energy  spectra  in  Figure  3.  7  show  a  broad  peak  at  2tt0  keV,  which 
is  the  approximate  position  of  th<*  most  probable  energy  loss.  The  av  rage 
energy  loss  is  about  3b0  keV  for  1,000  pm  of  Si  and  1-2  MeV  electrons.  The 
difference  is  due  to  the  large  amount  of  straggling  in  the  electron  energy  loss 
and  is  discussed  in  more  detail  in  Section  3.  3. 

Electron  anti-scattering  baffles  were  installed  as  shown  in  Figu-e  3.  , 

and  position  spectra  taken  with  PSD  1.  The  omplete  position  spectra  .re 
shown  in  Figure  3.8.  Comparison  with  Figure  3.7  shows  that  the  baffles  do 
not  affect  the  counts  in  the  region  of  beta  source  electrons,  except  for  the 
light  dip  in  the  Ru-Rh-106  Position  1  spectrum,  which  corresponds  to  a  baf¬ 
fle  location.  The  higher  energy  background  is,  however,  reduced  by  a  factor 
of  2  to  3  because  of  the  reduction  of  electron  scattering  in  the  magnet  gap. 

There  is,  however,  still  an  excess  of  electrons  in  the  1-2  MeV  range  just  a- 
bove  the  beta  source  end-point  energy. 

The  200  pm  detector  v  •.»  en  installed  to  achieve  the  full  configuration 
of  Figure  t.  1.  The  spectra  of  the  200  pm  detector  with  the  beta  sources  and  a 
50  keV  threshold  are  shown  in  f  igure  3.9.  The  peak  occurs  at  55  keV  and  is  » 
partial  artifact  of  th<  keV  threshold.  The  minimum  ionizing  energy  loss 
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Figure  i.  9  Spectra  of  the  200  pm  detector  with  the  Beta  Source# 
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for  1-2  MeV  electrons  is  about  70  keV,  so  the  most  probable  energy  loss 
(*55  keV  or  slightly  less)  is  significantly  lower  than  the  average  energy  loss. 
This  is  in  agreement  with  calculations  and  the  results  found  for  the  1,000  pm 
thick  PSD  1  and  is  discussed  in  more  detail  in  Section  3.  3. 

Spectra  of  PSD  1  in  coincidence  with  the  200pm  detector  with  a  50  keV 
threshold  gave  extremely  low  counts,  yielding  a  reduction  in  detection  effi¬ 
ciency  by  more  than  two  orders  of  magnitude.  To  achieve  good  detection 
efficiency,  it  was  necessary  to  lower  the  200  pm  detector  threshold  to  15  keV . 
This  was  partly  into  the  preamplifier  noise,  but  still  yielded  a  significant  im¬ 
provement  in  performance.  The  results  for  free  and  coincidence  operation 
of  the  PSD  1  are  shown  in  Figure  3.  10  for  the  Ru-Rh-lOb  beta  source  and  in 
Figure  3.  11  for  the  Sr-Y-90  beta  source.  These  tests  were  made  with  the 
source  and  two  detectors  in  a  straight  line,  with  a  small  collimator  in  front 
ot  the  detectors.  The  results  show  that  coincidence  operation  eliminates  most 
of  the  high  energy  pulses  and  the  uniform  background  in  position.  For  the 
straight  line  set-up,  the  detection  efficiency  is  reduced  by  about  a  factor  of 
10,  but  the  energy  and  position  spectra  narrow  significantly. 

The  free  and  coincidence  measurements  in  the  magnetic  electron  spec¬ 
trometer  configuration  are  shown  in  Figures  3.  12  and  3.  13  for  the  two  beta 
sources.  Comparison  with  the  Position  l  results  in  Figures  3.6  and  3.8  show 
that  the  addition  of  the  200  pm  detector  reduces  the  detection  efficiency  for 
1.6-2  MeV  electrons  by  about  a  factor  of  4  and  coincidence  operation  by  a- 
nother  factor  of  3.  The  overall  reduction  is  thus  about  a  factor  of  10.  Coin¬ 
cidence  operation  eliminates  most  of  the  spurious  counts  above  the  beta  source 
end* point  energy.  This  is  particularly  shown  by  the  strong  reductions  at 
3.  54  MeV  in  Figure  3.-12  and  at  2.  27  MeV  in  Figure  3.  13. 

The  coincidence  results  indicate  that  much  of  the  background  near  the 
beta  source  end  point  is  due  to  bremsstrahlung  generated  in  the  aluminum 
electron  collimator.  Since  coincidence  requirement  with  the  200  pm  detector 
eliminates  this  background,  most  of  it  is  generated  by  electrons  in  the  entrance 
part  of  the  collimator  in  front  of  the  200  pm  detector.  The  coincidence  re¬ 
quirement  thus  eliminates  most  of  the  background,  but  because  of  scattering 
in  the  coincidence  detector,  the  detection  efficiency  is  reduced  by  about  a 
factor  of  10. 


3.3  Analysis  of  Spectrometer  Response 


The  geometry  for  the  spectrometer  response  calculation  is  shown  in 
Figure  3.  14.  Since  the  geometric  factor  calculation  is  complex  and  involves 
some  uncertainty,  an  actual  flight  spectrometer  must  be  calibrated  with  elec¬ 
tron  beams.  For  the  present  case,  an  approximate  calculation  is  sufficient 
since  it  is  only  necessary  to  see  the  effect  of  various  parts  of  the  spectrometer 
on  the  geometric  factor. 


The  basic  geometric  factor  is  defined  by  A^,  A^ 

Go  .  A^/d2 


and  d,  and  is  given  by 
(3.  1) 
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The  input  vullimatur  nyal cm  admits  a  beam  of  total  angular  width  given  ap¬ 
proximately  by 


(3.2) 


From  the  collimator  exit  aperture  the  PSD  height  h  subtends  an  angle  of  ap¬ 
proximately 


,a"''  br) 


h  \  Tf  /  (3.3) 

Since  the  height  h  is  about  the  same  as  /A  , ,  only  about  hall  of  the  incoming 
electrons  will  actually  strike  the  sensitive  part  of  the  PSD  even  for  large  values 
of  .  In  the  r  plane,  there  is  no  electron  loss  because  of  focusing  by  the 

ISO  bend  in  the  magnetic  field.  Thus,  the  collection  factor  for  electrons 
leaving  A,  can  be  written  approximately  as 


F  (E)  =  0.  5  ,  6  >  tf 
o  ho 

=  0.5  0/  0  ,  <0  <  0 

h  o  h  o 


(3.  4) 


where  the  dependence  on  the  electron  energy  E  is  via  r  in  6^  in  (3.  3).  For 
a  strip  of  width  Ax  at  position  x  2r  ,  and  for  the  electrons  of  energy  E, 
the  count  rate  is  then  approximately 


AC  = 


G  F  (E)  Ax 

o  o _ 

2  (dr  /dE) 
c 


dJ  (E) 
c 


(3.5) 


where  d  J  (E)/d  E  is  the  incident  electron  flux  in  el/ (cm  -sr-sec-McV).  For 
the  test  magnet  assembly  and  electrons  above  about  1  MeV,  it  is  found  that 


_ 1 _  a  0.60  MeV/ cm  (3.6) 

2  (d  r  /  dE) 
c 

Using  the  values  oL,A^,  etc.  ,  for  the  test  magnet  assembly,  the  resulting 
values  of  GAE  (cm  -sr-McV)  for  Ax  =  1  cm  are  listed  in  Table  3.  1  for  1 

to  10  MeV  electrons.  The  values  are  for  AE  0.60  MeV,  as  obtained  from 
(3.  6). 

The  beta  source  tests  were  done  with  the  sources  about  8  1/2  inches 
(21.59  cm)  from  the  A,  aperture.  The  Sr-Y-90  source  was  measured  to  be 
10  mCi,  with  a  3/8  inch  diameter  aperture.  Taking  account  of  the  difference 
in  geometry,  this  yields  an  effective  source  strength  of  1.  12  x  10  el/  (cm  - 
sec-sr)  to  be  used  with  the  GAE  of  Table  3.  1.  For  the  Rn-Hh-106  beta  source, 
the  activity  was  measured  to  be  2.  5  mCi,  with  a  source  ^iameter  of  3/  16  inch, 
giving  an  effective  source  strength  of  2.80  >  10  el/(cm  -sec-sr). 

Source  measurements  made  in  the  air  require  an  additional  consideration 
for  electron  multiple  scattering.  Electron  multiple  scattering  and  the  resulting 
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Geometric  Factors  tor  Test  Magnet  Assembly 
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angular  and  radial  distributions  are  treated  in  Ref.  3.  1,  pp.  <18*1  -  29’.  lor 
electrons  penetrating  a  thickness  Mg/cni")  of  material  '/,  A(atomic  number 
and  weight),  the  probability  distribution  in  angle  is 
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with  v  the  electron  velocity,  c  the  velocity  of  light  and  p  the  electron  momen¬ 
tum  with  (pv)  in  MeV.  The  lateral  displacement  distribution  after  traveling  y 
cm,  which  contains  the  t  g/  cm  of  material,  is 


F(r,  y)  dr 


brdr 


*3r*”  /  (<W*’>  y“) 


<0  >  y 


where<W  >is  calculated  for  t  with  (3.8). 


(3. 


The  above  distributions  must  be  used  twice,  once  for  the  path  from 
the  beta  source  to  A,,  and  once  from  A,  to  the  PSD,  ,Kor  the  first  path, 
the  parameters  art;  y  -  21.59  cm,  t  =  0.0259  g/  cm  ,  and  R  0.318  cm. 
The  results  for  “  >  and  P(r  <  R)  a  re  given  in  Table  3.2,  where  P(r  <  R)  is 
obtained  from  integration  of  (3.9)  and  is 

P(r  ^  R)  -  1  -  e‘3R  /  ^  ^  V  ] 


(3.  10) 


For  the  second  path,  the  calculation  is  somewhat  more  complex  because  to 
first  order  the  spread  in  tin  r  plane  can  be  neglected,  and  only  those  elec¬ 
trons  scattered  out  of  this  plane  and  beyond  _+  h/ 2  are  lost.  This  can  be  ap¬ 
proximated  by  taking  the  square  root  of  (3.  10),  splitting  the  circular  distri¬ 
bution  in  half.  Doing  this  gives  probabilities  ranging  from  0.97  at  1.5  Me V 
to  0.77  at  5.0  MeV,  and  they  arc  given  in  the  last  column  of  Table  3.2. 

The  measured  position  spectra  are  converted  to  incident  tlux  by 


d  J  (E) 

o 

d  E 


_ (Measured  counts/ channel)  _ _ 

(cm/ eh)  UTKCAE)  P(rS0.  '18)1  (hi 


(.11) 


where  GAF  is  taken  from  Table  3.  1,  P  and  P  from  Table  3.2,  and  (cm/ -h) 
and  (AT)  (count  time  in  seconds)  are  tor  the  MCA  spectrum  in  question.  I  he 
spectra  In  Figure  J.  6  have  (cm  ch)  0.0 (  '  and  (AT)  B  x  60  I  nd  .  witl 
the  result!  tret  .  tra  being  git  n  in  Table  3. 3.  The  measured 

spectra  are  compared  t"  the  calculated  spectra,  based  on  lh<  Sr-\-90  and 
Ru-Rh-106  spectra  given  in  Ref.  3.2.  Consid*  the  approximations  and 

corrections  made,  the  agreement  is  quite  good.  I  he  dip  in  (measure  d/ cal¬ 
culated)  for  Ru-Rh-106  at  2.75  MeV  may  be  duo  to  the  dip  in  magnetic  field 
at  the  gap  between  the  first  and  second  row  of  magnets  (see  i  igures  2.5  and 
2 .  b) ,  since  the  circular  path  tor  electrons  has  a  maximum  length  in  this  dip 

for  about  2.  75  M«V  (••«  Figure  2.  7).  (A  tin  »r  *lf<  tv  ild  e  expect*  i  ■  r 

b  -  b.  5  MeV  electrons,  but  no  electrons  ..{  this  ene  rgy  art  emitte  d  b\  the  a  ta 
sources  used.)  The  rise  in  (mea  surrd/  cal*  ulate-d)  at  the*  beta  source  end¬ 
point  energy  and  the  existence  of  counts  beyond  the  end  point  are  due  tei  brern- 
sstrahlung  background  and  possibly  some  contribution  from  the  scattered  rl<  .  - 
irons.  The  measured  spectra  in  f  igure  3.  b  thus  show  good  agreement  with 
that  calculated  for  the  beta  sources  but  show  an  important  background  at  higher 
e*  no  r  g  i  e  s . 

The  addition  of  the  anti-scattering  battles  has  only  a  minor  effect  on  the 
beta  spectra  but  does  reduce  the  higher  e*nergy  background  by  a  fai  tor  of  2-3, 
as  shown  in  Figure  3.8.  The  dip  in  the  Ku-Rh-106  spectrum  between  2  and 
MeV  is  accentuated  somewhat  because  now  a  baffle,  which  does  eliminate  some 
non-scattered  electrons,  adds  to  the  magnetic  field  dip  effect  at  2.  75  MeV. 

It  is  expected  that  a  flight  unit  magnet  assembly,  which  would  not  have  the 
gap  between  the  magnet  row  s  (see  Section  4),  would  have  a  more  uniform  mag¬ 
netic  field  and  thus  not  show  significant  non-uniformities  in  spectral  response. 
Any  non-uniformities  that  remain  will,  however,  not  be  important  since  the 
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TABLE  3.2 


Electron  Transmission  after  Multiple  Scattering  In 

Air 


Electron 
Ene  rg  y 
(Me  V) 

(v/c)“ 

(pv)“ 

(MeV) 

</>a 

(radians)" 

P(r<0.  318)b 

Pj(h)C 

=  v^P^r  <  0.  35) 

1.5 

0. 9354 

3.  539 

0.02687 

0.0239 

0.  970 

2.0 

0. 9586 

5.  794 

0.  01634 

0.  0390 

0.  939 

3.0 

0. 9788 

1 1. 810 

0.00799 

0. 0782 

0.  888 

5.0 

0. 9914 

29.  85 

0.00315 

0. 1867 

0.  770 

a  ,  2 

Calculated  for  0.  0259  g/  cm  of  air  (Z  7.22,  A  14.  485)(Eq.  (3.  8)) 
b  ,  2V 

Calculated  for<6  )  and  y  21.  59  (Eq.  3.  10). 

Calculated  tor  a  path  length  nl  »r  ,  with<*^  >  va rying  with  this  path  length. 
_  e 

See  text. 

TABLE  3.  3 


Comparison  of  Measured  and  Calculated  Beta  Source 

Sped  ra 


Electron 

dJ  (E)/dE 
c 

in  el/  (cm*"  - 

sec  -  sr  -  MeV) 

Ene  rg  y 

(MeV) 

Calc  vilated 
IRcf.  3.2) 

Measured 
(P»U-  )-  bj 

/  Measured  \ 
\  Calculated  / 

Calculated 
LH.  t,  3,  2J 

Measured 

(E»Ji.  b) 

/  Nlea  sured 

\Calc  ulated 

1. 75 

3.  47  +6 

2.41  4  6 

0.  70 

1. 20  4  6 

1.  27  4  6 

1.06 

2.  00 

1.23  4  6 

1.51  4  6 

1. 23 

1.01  4  6 

1.01  4  6 

1.00 

2.25 

0.  00 

4.  94  4  5 

.... 

8.  12  4  5 

6.  40  4  5 

0.  79 

2.  50 

— 

1.66  4  5 

.... 

5.  88  4  5 

4.  33  4  5 

0.  74 

2.75 

3.  92  4  5 

2.  36  4  5 

0.60 

3.00 

1.  99  4  5 

1.69  4  5 

0.  84 

3.25 

5.  60  4  5 

1.1445 

2.04 

3.  50 

0.00 

7.  07  4  4 

.... 

*  3.  47  4  6 

=  3.  47  x 

106,  etc. 
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flight  unit  would  bo  calibrated  directly  with  electrons  from  an  accelerator. 

The  electron  energy  loss  spectra  in  the  1,000pm  thick  PSD  are  shown 
in  Figure  3.7  and  show  a  peak  near  275  feV.  For  1.6  to  4.4  MeV,  using  the 
stopping  powers  for  Si  in  Ret.  3.2,  electrons  should  deposit  an  average  ener¬ 
gy  of  3o0  keV.  From  Ref.  3.  1,  pp.  252  -  2b 7,  it  is  found  that  for  relativistic 
electrons  the  most  probably  ei\e rgy  loss  can  be  written  as 

/  mc"  P  *  ,  \ 

AE  =  W  In  - 7 - ;  ]  -  P"  -  K  (3.12) 

p  1  V  ru+r)  / 


who  re 


0.  3006  /.  t 

A 


(3.  13) 


with  me  the  electron  rest  mess  in  r  \  ,  /?  v/ c  is  the  electron  velocity  di¬ 
vided  by  the  velocity  of  light,  1  is  thy  average  excitation  potential  (in  eV)  of 


material  /,  A,  of  thickness  t(g/cm"),  mc^K*  0.  355  is  a  constant.  For  1,000pm 
of  Si  (X  14,  A  28.086,  t  0.  233  g/cm  )  and  2  MeV  electrons  this  gives 


2iE  (2  MeV  el)  2  70  keV 
P 


(3.  14) 


which  is  in  close  agreement  with  the  measurements.  As  already  mentioned 
when  disc  us  sine  figures  3.4  and  3.5  (Section  3.2),  the  breins  st  rahlung  back¬ 
ground  appears  to  have'  a  higher  average  energy  than  the  direct  electrons,  so 
a  100  -  150  keV  threshold,  well  below  the  275  keV  peak,  gives  the  best  signal/ 
noise  ratio.  A  threshold  below  100  keV  will  not  further  improve  this  ratio. 

Use  of  a  t)unn<  r  PSD  would,  however,  require  a  lower  threshold  and  introduce 
more  noise  in  the  position  signal  calculation.  The  measured  electron  energy 
loss  distributions  are  also  in  reasonable  agreement  with  the  calculations  of 


Ref.  3.4. 


The  final  configuration  uses  a  200pm  Si  detector  in  the  entrance  colli¬ 
mator  to  provide  a  coincidence  signal.  The  average  energy  loss  for  1.6  to 
4.4  MeV  electrons  is  about  71  keV  (from  data  in  Ref.  3.3),  while  the  most 
probable  energy  loss  from  (3.  12)  is  48  keV.  Energy  loss  measurements  with 
a  50  keV  threshold  show  the  most  probable  energy  loss  to  be  about  50  keV  or 
less  (Figure  3.9)  and  are  thus  in  agreement  with  calculations.  Actual  coin¬ 
cidence  measurements  with  the  200  pm  and  PSD  detectors  showed  that  a  15  keV 
threshold  on  the  200  pm  detector  is  necessary  to  give  a  reasonable  detection 
efficiency.  For  the  straight  line  arrangement,  the  maximum  coincidence 
detection  efficiency  is  only  about  lO^c  and  is  due  to  the  multiple  scattering  of 
electrons  in  the  200  pm  detector.  The  multiple  scattering  in  the  200  pm  de¬ 
tector  spreads  the  electron  beam  so  that  only  about  lOTt  of  those  producing 
large  pulses  pass  through  the  collimator  in  front  of  the  PSD. 

Use  of  the  200pm  detector  in  the  magnet  assembly,  with  the  complete 
set-up  as  shown  in  f  igure  3.  1,  gives  the  results  shown  in  f  igures  3.  12  and 
3.  13.  A  rough  calculation  of  the  multiple  scattering  effect  ol  the  200pm  de¬ 
tector  gives  a  reduction  of  about  a  factor  of  10.  The  measured  reduction  is 
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about  a  factor  of  4  for  no  coincidence  and  an  additional  factor  of  3  with  coin¬ 
cidence.  The  factor  of  4  is  considered  reasonable  agreement  with  the  cal¬ 
culated  10,  since  the  complete  scattering  effect  is  difficult  to  calculate  be¬ 
cause  the  200  pm  detector  scattering  adds  on  to  the  air  scattering.  The  tactor 
of  3  coincidence  reduction  is  less  than  the  factor  ol  10  for  the  straight  line 
geometry.  The  coincidence  reduction  results  from  the  correlation  between 
angular  scattering  and  energy  loss.  Electrons  lose  energy  by  scattering, 
so  those  electrons  which  go  through  the  200  pm  detector  on  the  straightest 
paths,  and  hence,  are  most  likely  to  reach  the  PSD,  also  are  likely  to  have 
the  lowest  energy  losses  in  the  200  pm  detector. 

The  net  result  of  adding  the  coincidence  detector  is  to  reduce  the  geo¬ 
metric  factors  by  nearly  a  factor  of  10.  It  also  reduces  the  higher  energy 
background  by  more  than  order  ot  magnitude  and  so  improves  the  signal/ 
noise  ratio.  The  spectra  in  Figures  3.  12  and  3.  13  show  significantly  im¬ 
proved  behavior  at  the  beta  end-point  energies. 

A  minor  problem  occurs  with  the  low  energy  threshold  (1  krV)  re¬ 
quired  for  the  200  pm  detector.  This  is  about  the  limit  ol  what  can  be  »chiev  d 
with  present  room-temperature  preamplifiers.  Use  of  a  thicker  detector 
would  raite  this  threshold,  but  it  w  .>uld  also  increase  multiple  scattering 
effects;  therefore,  no  increase  in  detection  efficiency  i*  111  lv  to  result.  The 
200  pm  detector  is  quite  probably  near  the  optimum  thickne  s. 

Placing  the  200  pm  detector  nearer  to  the  collimator  exit  aperture  (A, 
in  Figure  3.  14)  would  increase'  th<  efficiency  somewhat  since  the  magnetic 
focusing  in  the  r  plane  ■a  aid  r.  .liify  ne  dirr  tion  of  scatl<  ring.  This 
would  also  increase  the  number  of  electrons  striking  the  magnet  pi  lea  and 
so  increase  background  from  b  r>*m  s  s  t  rahlung  and  scattered  eli  tr  ns.  This 
would  thus  not  be  likely  to  increase  the  signal/ noise  ratio  for  detecting  high 
energy  electrons  in  the  presence  of  low  energy  electrons. 

The  final  configuration.  Figure  3.  14,  is  thus  close  to  the  best  conligur- 
ation  for  the  Magnetic  Electron  Spectrometer.  With  additional  lead  shielding 
for  reduction  ot  b  reins  st  rahlung  background,  and  a  possible  scintillator  be¬ 
hind  the  PSD'  S,  this  is  the  final  design  configuration,  on  which  the  prelimi¬ 
nary  flight  unit  design  is  based. 

4.  PRELIMINARY  DESIGN  OF  A  FLIGHT  UNIT  MAGNETIC  ELECTRON 
SPECTROMETER 

The  results  in  Section  3  allow  a  preliminary  design  for  a  flight  unit  mag¬ 
netic  electron  spectrometer  to  be  made.  The  properties  of  this  flight  unit  can 
be  specified  with  reasonable  accuracy.  The  flight  unit  magnet  arrangement  is 
shown  in  Figure  4.  1.  The  design  uses  the  same  iron  yoke  design  as  the  test 
magnet  assembly  but  uses  some  extra  magnets  to  improve  the  measurement  at 
10  MeV.  The  design  uses  19  -  1  1/  1  b  x  1  1  /  1 6  x  1/ 2  inch,  3  -  1 1/  16  x  1/2  x 
1/2,  and  2  -  1/2  x  1/2  x  l/2  inch  magnets  on  each  plate.  The  magnets  are 
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Figure  4.  1  Magnet  Arrangement  for  Flight  Unit  Magnet  Assembly. 


selected  from  a  larger  number  on  the  basis  of  magnet  field  measurements 
to  produce  a  more  uniform  field.  The  largest  magnets  have  a  tab  ground 
off  on  one  corner  to  allow  the  use  of  holding  screws  during  and  after  assem¬ 
bly.  Elimination  of  the  screw  gaps  (Figure  2.2)  should,  in  combination  with 
magnet  selection,  allow  the  magnet  assembly  to  have  an  average  field  in  the 
gap  of  4.  50  _+  0.  15  kg,  a  _+37o  uniformity. 

The  paths  of  1  to  10  MeV  electrons  in  the  magnet  gap  are  shown  in  Fig¬ 
ure  4.2.  There  is  sufficient  edge  field  to  allow  magnetic  focusing  even  for 
the  10  MeV  electrons.  With  the  field  uniform  to  _+  3%,  deviations  from  the 
paths  in  Figure  4.2  should  not  be  significant  at  the  performance  level.  The 
effects  of  fringing  fields  will  cause  slight  changes  in  the  curved  and  straight 
paths  at  the  magnet  pole  edges  but  should  have  only  a  small  effect  on  the  exit 
positions  for  1  to  10  MeV  electrons. 

The  complete  entrance  collimator,  detector,  and  baffle  design  is  shown  in 
Figure  4.  3.  The  entrance  collimator  contains  a  front  magnet  assembly  to 
shield  the  coincidence  detector  from  electrons  below  100  -  200  keV  and  so 
avoid  pile-up  problems  from  large  fluxes  of  low  energy  electrons.  The 
entrance  collimator  is  made  from  aluminum  to  reduce  bremsstrahlung  prod¬ 
uction.  The  200  Hm  surface  barrier  detector  has  a  thin  aluminum  foil  light 
shield,  which  adds  negligibly  to  the  electron  scattering  of  the  detector  itself. 

The  detector  at  the  magnet  exit  line  is  designed  to  be  a  strip  surface 
barrier  detector,  1 , 000  gum  thick,  and  segmented  into  12  individual  rectan¬ 
gular  detectors.  The  detector  strip  need  not  be  a  single  piece  of  silicon  but 
must  be  mounted  with  only  a  small  dead  section  between  adjacent  strips.  The 
detector  array  is  shielded  from  behind,  top,  bottom,  and  from  the  entrance 
collimator  by  lead  to  reduce  the  brems strahlung  background.  A  shielded  de¬ 
tector  at  the  high  energy  end  provides  a  background  measurement. 

The  strip  detector  is  designed  with  7  segments  of  1  cm  width  at  the  low 
energy  end,  4  segments  of  2  cm  width  at  the  high  energy  end,  and  a  1  cm  wide 
segment  for  background  measurement.  The  junctions  between  adjacent  seg¬ 
ments  have  a  0.25  inch  long  strip  of  1/  16  inch  copper  extending  to  the  magnet 
pole  edges  as  an  electron  shield.  Five  baffles  of  30  degree  arc  length  are  in¬ 
stalled  as  shown  in  Figure  4.  3.  The  baffles  are  limited  to  30*  to  reduce  their 
interference  with  the  magnet  focusing  properties. 

The  segmented  strip  detector  design  was  chosen  over  the  use  of  position 
sensitive  detectors  for  a  number  of  reasons.  The  use  of  position  sensitive 
detectors  would  have  required  at  least  three  (3),  which  requires  six  (6)  pre¬ 
amplifier/ amplific  r  chains,  already  half  of  the  requirement  for  the  segmented 
strip.  The  position  calculation  electronics  are  easily  the  equivalent  of  the  re¬ 
maining  six  (6)  preamplifie r/  amplifier  chains;  therefore,  the  segmented  strip 
does  not  require  a  significant  amount  of  additional  electronic  components. 

The  segmented  strip  detectors  can  be  made  to  operate  with  a  total  re¬ 
sponse  time  of  less  than  one  microsecond,  while  the  position  calculation  elec¬ 
tronics  typically  requires  several  microseconds.  Thus,  the  segmented  strip 
can  handle  count  rates  (and  fluxes)  nearly  an  order  of  magnitude  higher  than 
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in  Flight  Unjt  Magnet  Assembly 


baffles  (5  total) 


the  position  sensitive  detectors. 

Finally,  the  segmented  strip  can  operate  with  all  segments  in  anti- 
coincidence,  thus  strongly  reducing  the  effect  of  backscatte red  electrons 
from  a  low  to  higher  energy  position.  A  position  sensitive  detector  cannot 
eliminate  such  backscatte  ring  events  on  its  own  length.  The  segmented 
strip  should  thus  provide  much  sharper  cut-offs  for  electron  spectra  vary¬ 
ing  strongly  with  energy. 

The  position  sensitive  detectors  can  in  principle  give  better  electron 
energy  resolution  via  the  position  signal.  This  is  somewhat  deg  raded  by  the 
electronic  noise  added  to  low  electron  energy  losses  and  the  scattering  of 
electrons  in  the  silicon.  Ultimately,  it  is  limited  by  the  number  of  energy 
channels,  which  is  determined  by  the  number  of'scalers  used.  The  design 
in  Figure  4.  3  uses  12  scalers  for  an  1  1  point  energy  spectrum  and  is  un¬ 
likely  to  be  improved  significantly  without  a  large  increase  in  electronic 
components . 

The  basic  electronics  for  the  electron  spectrometer  is  outlined  in  figure 
4.4.  Counting  the  200  ►cm  detector,  a  total  of  13  preamplifier/ amplifie  r  chains 
are  required.  The  basic  logic  for  coincidence  (200  pm  detector)/ anti -coinci - 
dence  operation  is  shown.  A  set  of  12  scalers  count  the  II  electron  channels 
and  one  background  channel.  Additional  circuitry  for  control  of  cycle  times, 
etc.,  are  not  included  as  they  would  depend  on  the  interface  requirements 
for  the  satellite  in  which  the  electron  spectrometer  is  placed. 

Figure  4.3  also  shows  a  space  for  an  anti -coincidence  scintillator  be¬ 
hind  the  segmented  strip  detector.  This  would  allow  the  use  of  the  spectro¬ 
meter  without  the  200  pm  detector  and  in  the  presence  of  penetrating  radiation. 
Because  of  the  large  geometric  factor  of  the  scintillator,  the  penetrating  flux 
(generally  high  energy  protons)  cannot  be  too  great  before  the  scintillator  count 
rate  becomes  too  high  to  be  usable.  The  anti -coincidence  scintillator  is  thus 
indicated  only  as  a  possibility  for  certain  very  specific  flux  measurement  con¬ 
ditions. 

The  properties  of  the  magnetic  electron  spectrometer  electron  channels 
are  listed  in  Table  4.  1.  The  GAE  factors  are  listed  both  with  and  without  the 
200  Hun  coincidence  detector  in  place,  the  difference  being  the  electron  scat¬ 
tering  factors.  The  GAE  values  in  Table  4.  1  are  estimated  to  be  accurate  to 
a  tactor  of  3.  In  practice,  a  flight  unit  spectrometer  would  lie  calibrated 
with  electron  beams  to  reduce  the  uncertainties.  This  would  also  calibrate 
the  energy  response  of  each  channel.  The  GAE  values  for  no  200  pm  detector 
hold  for  operation  with  or  without  the  anti -coincidenc e  scintillator. 

The  physical  properties  of  the  preliminary  flight  unit  design  are  given 
in  Table  4.2.  The  magnet  assembly  weight  includes  all  iron  pieces,  the 
magnets,  and  the  aluminum  braces  (see  Figure  2.3).  The  listed  design  in¬ 
cludes  the  preamplifier/ amplifier  chains  with  the  magnet/detector  assembly 
and  a  separate  electronics  unit  for  power/  signal  processing.  If  the  electronics 
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TABLE  4.  1 


Electron  Detection  Properties  of  the  Flight  Unit  Design 


Electron  Electron  Average  Bin  GaE  with  GAE  without 


Bin 

Numbe  r 

Energy  Range 
(Me  V ) 

Energy 

E  (MeV) 

Width 

AE  (MeV) 

200  un  detector 
(cm  -  s  r  -  MeV) 

200  tym 
(cm  - 

detector 
s  r  -  MeV) 

l 

0.87  - 

1.51 

1.  19 

0.  64 

7.  1 

- 

6 

5.  5 

- 

4 

2 

1.51  - 

2.  17 

1. 84 

0.  6b 

1.5 

- 

5 

5.  3 

- 

4 

3 

2.  17  - 

2.  84 

2.  51 

0.  67 

2.  2 

- 

5 

4.  5 

- 

4 

4 

2.84  - 

3.  51 

3.  18 

0.  67 

2.  8 

- 

5 

3.  9 

- 

4 

5 

3.51  - 

4.  18 

3.  85 

0.  07 

3.  3 

- 

5 

3.  4 

- 

4 

6 

4.  18  - 

4.  85 

4.  52 

0.  67 

3.  8 

- 

5 

3.0 

- 

4 

7 

4.  85  - 

5.  52 

5.  19 

0.  67 

4.  1 

- 

5 

2.  7 

- 

4 

8 

5.  52  - 

6.  87 

6.  20 

1.  35 

8.  5 

- 

5 

4.  8 

- 

4 

9 

6.  87  - 

8.  22 

7.  55 

1.  35 

8.  7 

- 

5 

4.  2 

- 

4 

10 

8.22  - 

9.  57 

8.  90 

1.  35 

8.  9 

- 

5 

3.  5 

- 

4 

1 1 

9.  57  - 

10.  92 

10.  25 

1.  35 

9.  1 

- 

5 

2.  9 

- 

4 

12  Background 


*  7.  1  -  6  =  7.  1  x  10‘6,  etc. 

TABLE  4.  2 


Estimated  Physical  Properties  of  Flight  Unit  Magnetic 
Electron  Spectrometer 

Si/e:  Magnet/ Detector  Assembly  -  8x9x5  inches 

Electronics  Unit  -  6  x  6  x  10  inches 


Weight:  Magnet/ Detector  Assembly 

Magnets,  iron  yoke,  braces 
Collimator  ♦  lead  shielding 

Detectors  +  preamphfie  r/ amplitie  r  chains 
Electronics  Unit 

TOTAL  = 


-  19  lbs. 

-  3  lbs. 

-  2  lbs. 

-  6  lbs. 

30  lbs. 


Power  Consumption  -  from  spacecraft  buss 


10  W 


Note:  If  the  electronics  unit  is  combined  with  the  magnet/detector  assembly, 
the  combined  size  would  be  about  8x9*  10  inches,  and  a  small  weight  saving 
would  result.  The  actual  design  of  2  or  1  housings  can  be  done  to  fit  space¬ 
craft  requirements. 


42 


urut  can  be  combined  with  the  magnet/ detector  assembly,  a  small  savings 
in  weight  could  be  achieved.  This  latter  design  would  be  preferable  from 
an  electronics  design  viewpoint  as  it  eliminates  external  signal  cables  for 
the  detectors.  The  actual  use  of  2  or  1  housings  can  be  done  according  to 
spacecraft  requirements.  The  power  consumption  of  10W  is  input  power 
to  the  DC/DC  converter  from  the  spacecraft  buss  and  covers  all  instru¬ 
ment  requirements. 

5.  CONCLUSIONS 

A  preliminary  design  of  a  flight  unit  magnetic  electron  spectrometer 
for  less  than  1  to  greater  than  10  MeV  electrons  has  been  described.  The 
design  is  based  on  an  extensive  series  of  tests  with  a  rare-earth  cobalt  mag¬ 
net  assembly  and  position  sensitiv<  solid  state  detectors.  These  tests  proved 
that  a  magnet  assembly  could  be  built  to  analyze  to  at  least  !0  Me  V  electrons. 
The  major  background  sources  were  found  to  be  internally  scattered  electron:  , 
and  brents st rahlung ,  both  internail)  and  externally  generated.  For  flight  in¬ 
to  the  inner  belt,  background  du<  to  penetrating  protons  woulu  also  hav<  to  b« 
conside  red. 

The  final  configuration  uses  a  lead- shielded  segmented  - s i  r: p  surlace 
barrier  detector  at  the  magnet  exit  plane  and  a  200  pm  Iran  mission  detector 
in  the  entrance  collimator  for  coincidence  measurement.  Scattered  electron 
detection  is  reduced  by  baffles  in  the  magnet  assembly.  I'he  segmented- 
strip  detector  arrangement  was  selected  because  it  is  no  more  complex  than 
the  position  sensitive  detector  a;  rangement  but  allows  highi  r  count  r.^tes  and 
provides  anti -coincidence  discrimination  against  low  energy  electrons  which 
are  backscattc  red  to  a  higher  energy  position. 

The  final  design  uses  cithe^  two  or  one  housings  with  a  total  volume 
of  about  8  x  9  x  10  =  720  inches  ,  a  total  weight  of  about  30  lbs.  ,  and  a  total 

power  consumption  of  10  W. 
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Appendix  A 

The  I'rajectory  of  a  Charged  Particle  in  the 
Fringing  Field  of  a  Magnet  for  Arbitrary  Incidence  Angle 


A.l  Introduction 

It  is  often  necessary  t-  find  t:  <  tra-ectory  of  a  particle  in  a  magnet's  fring 
ing  field,  similar  ti<  that  shown  in  Fig.  A.l.  Here 


h(y ) 


Ihil 


H 

o 


(A.l  ) 


where  H(y)  is  field  -  •  .  «*  .statu  \  (perpendicular  to  the  pole  face)  from 

a  point  inside  the  p  1<  'und.irv  a*  w  ■  n .  h  H(yl  may  be  assumed  to  have  its  uni¬ 
form-field  value  H  .  With  referent  .  *  Fig.  A.  I,  the  general  problem  is  that  of 

o 

finding,  for  known  h(y>,  the  trajectory  \<y)  and  the  angle  tty)  made  by  a  tangent 
to  the  trajectory  with  the  uniform  field  boundary. 

The  only  value  of  «|y>  that  is  normally  f  significance  is  f(y  0).  How¬ 

ever,  x(y)  is  often  o:  interest  for  arbitrary  y.  in  order  to  determine  whether 
the  trajectory  intercepts  collimators  etc.  located  at  specific  values  of  y. 

The  approach  taken  here  for  arbitrary  incidence  angle  •  .  is  similar  to  that 
of  Coggeshell  (P  ef.  A.l)  and  Coggeshall  and  Mu  skat  (R  ef.  A.  2 )  for  normal  inci¬ 
dence  (s.  0).  T  is  assumed  tha*  H  is  a  function  H(y);  it  does  not  depend  on  x. 
Clearly,  this  simplifies  the  problem  considerably.  In  a  realistic  case,  it  is  a 
good  approximation  if  the  distance  x  f  at  which  the  integration  must  be  started 
is  small  compared  to  the  x-dimension  of  the  magnet  pole  faces.  Of  course,  in 
order  to  place  the  results  in  an  external  reference  frame  it  is  necessary  to 
specify  y  .  the  uniform  field-pole  boundary  distance.  As  seen  in  Fig. A.l,  for 
that  particular  magnet  (pole  faces  small  compared  to  gap  width;  yoke  sufficient 
to  contain  external  field  quttr  well)  this  distance  was  about  1.  6  gap 
widths.  Hence,  in  addition  to  h(y),  it  will  be  assumed  that  Yp  is  also  specified. 
The  field  reversal  shown  in  Fig. A.l  is  caused  by  the  fact  (hat  the  entrance 
aperture  of  the  assembly  penetrates  the  yoke. 
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The  integration  is  started  at  the  point  xjc(yac),  where  the  particle  trajec¬ 
tory  makes  the  initial  angle  with  *he  y-axis.  From  simple  geometrical  con¬ 
siderations,  x .  (y )  can  be  found.  Then  Ax(y).  the  fringing  field  induced  devia¬ 
tion,  can  be  added  to  obtain  x(y) 

x(y)  =  x.(y)  +  Ax(y)  (A.2) 

=  x.  +  ytanc  +Jkx(y)  (A. 3) 

10  '  i 


w  h  e  r  e 


(A. 4) 


in  terms  of  the  chosen  initial  conditions  x  ,  v  .  and  .  ..  Thus,  in  terms  of 

*  *  i 

these  conditions  and  the  deviation  Ax(y).  the  trajectory  x(y)  can  be  found  from 
(A.  3). 

Conceptually,  it  is  convenient  to  consider  the  problem  as  above.  However, 
in  the  procedure  below  x(y)  is  determined  directly,  as  is  <(y).  A  particular 
analytical  function  for  h(y)  is  then  shown  to  fit  the  curve  in  Fig.A.l.  This  func¬ 
tion  is  used  in  the  general  equation  for  *(y),  allowing  an  entirely  analytical  re¬ 
sult  for  that  quantity.  It  is  found  that  this  then  yields  a  result  r«-quiring  only  one 
simple  numerical  integral  for  xly). 

A. 2  General  Equations  for  x(y)  and  *(y). 

The  differential  equations  for  a  charged  particle  of  initial  velocity  V(  and 

radius  of  curvature  r  in  the  median  plane  uniform  field  (inside  magnet)  are 

o 

-v 

x  =  ~  h(y)y  (A. 5) 

o 

V 

y  =  —  h(y )  x  (A. 6 ) 

o 


Multiplying  the  first  by  x, 


the  second  by  y.  yields 


or 


-xx  =  yy 


dt  dt 
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Thu# 


2 


(A. 7) 


.  2  .  .  2  _ 
x  +  y  =  v 


and  (choosing  the  proper  sign  for  y) 


V-*T 


Given  x.  y  i»  then  found  from(A.8).  From  (A. 5) 

y 


?f 


x(y )  -  —  /  h(y')dy'  +  x_ 

o 


(A. 8) 


(A. 9) 


Thus,  if  u.e  define  the  quantity 


R(y )  -  -  j  h(y  )dy '  =  J  h(y')dy' 


(A. 101 


1  hen 


x  —  R(y)  -  v  mni 

r  o  i 

o 


Division  <>f  (A.ll)by  (A. 8)  yieldi 


x  dx 

y 


r  sin» .  -  R (  y ) 
O  l 


Un,,v)  \r  2  -  (r  Sin,.  -  Riy))1!* 
o  Vo  t 


(A.  1  1  ) 


(A. 12) 


This  equation  thus  defines  *(y).  as  sought.  If  the  analytic  expression  chosen  to 
fit  h(y)  is  integrable  in  (A. 10),  #  |y)  is  defined  analytically  from  (A.  1  2).  Although 

tan<(y)  is  required  to  find  x(y)  below,  the  angle  cq  at  entry  to  the  uniform  field  is 
more  easily  determined  -  as  ran  be  seen  from  (A. 12)  -  from 


sin  *o  sin  «.  -  R(0)/ r^ 


(A. 13) 


x(y)  can  now  be  found  by  integration  of  (A. 12).  which  must  be  carried  out  num¬ 
erically  for  any  reasonably  accurate  h(y)  representation.  The  result  is 
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In  summary,  fur  any  h(y),  R(y)  is  defined  by  (A.  10).  «  and  x  are  found 

'  o  o 

from  i A.  1  3)  and  (A.  1 <>).  and  the  trajectory  x(y )  from  (A. 1 5)  w  ith  tan»(y )  obtained 
from  (A.  1  2). 

A.  3  Results  for  Rfy)  with  an  Analytical  Approximation  for  h(y). 

Let  it  be  assumed  that  h(y)  can  be  written 

N 

i 

h|y)  f'^^ay1'1  (A. 17) 

j=l  J 

That  thiB  is  a  good  approximation  for  the  dependence  of  h  on  y  in  Fig.  1  is  shown 
in  Section  A. 5  below.  The  number  of  terms  N.  that  must  be  used  is  *  9.  In 
order  for  hfy)  to  have  the  value  unity  at  y  0,  we  must  have  aJ  1.  A  further 
condition  is  obtained  by  requiring  that  dh/dy  0  at  y  0,  which  requires  K'a^. 

These  conditions  are  used  in  Section  A. 4  in  development  of  a  least  square  pro¬ 
cedure  for  determining  the  constants  a  that  ’best  fit"  an  experimental  hfy) 
curve  such  as  Fig. A.  1. 
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Her e,  (A.  17)  is  substituted  into  (A.  10)  and  integrated.  After  some  re¬ 
arrangement  the  result  can  be  written  as  follows,  assuming  y^— 


(A. 18) 


N. 


i 


(A. 19) 


Thus,  given  all  of  the  a  and  K.  the  g  can  be  determined  from  (A.  19).  Insofar 
as  the  trajectory  determination  is  concerned,  further  reference  t"  the  a.  ne,*d 
not  subs  ,-quently  be  made,  since  only  the  value  of  Ry )  from  (A.  18)  i»  needed 
in  (A. 12)  for  use  in  (A.  1  5).  Note  that 

RIO)  -  g t / K  (A. 20) 


sc  that  *hr  entrance  angle  c  is  easily  found  from  /A.  13)  oner  the  h(y)  fit  (A. 17) 

o 

has  been  completed. 

A. 4  Least  Square  Procedure  for  Fitting  hty). 

Because  hly)  contains  a  non-linear  parameter  K,  the  usual  least  square 
procedure  must  be  modified  somewhat.  It  is  not  difficult  to  find  the  general 


range  of  values  that  K  must  have,  given  the  experimental  data  )v,  y^,  i  1, 
that  must  be  fitted.  Here  we  shall  have,  as  noted  above, 


(A. 21 ) 


There  are,  therefore,  N  -  1  coefficients  to  be  determined  from  the  N  values 

»  P 

of  h(N  >  N  ). 

P  * 
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In  the  case  of  interest  here,  then, 


h(y )  -  e 


a  y"*  1 ;  a  =  1 

y  1 


(A .22) 


Because  the  deflection  depends  on  the  man  nit  tide  of  h(y),  the  trajectory  that 
deviates  least  from  that  which  would  actually  occur  in  the  field  h(y)  should  be 
obtained  by  minimizing  the  actual  differences  in  field  magnitudes,  rather  than 
minimizing  the  fractional  differences.  Thus,  we  write  the  sum  of  the  squared 


deviations  as 


N 

£ 


h(y. )  -  h. 


(A. 23) 


The  procedure  used  here  is  to  vary  a  throughout  a  selected  range  of 

values.  For  each  given  a,  the  values  of  a  (j  3,  N.)  are  determined  that  pro- 

1  .  1  1  2 

vide  the  best  least  square  fit  to  (A. 23),  and  the  associated  value  A(a  )  is  found. 

1  2 

The  optimum  value  of  a,,  and  associated  a  (j  -  3,  N.),  is  that  for  which  A(a,) 

2  j  '  \  2 

is  minimum. 

For  convenience  we  let  c  ,  a  ,  i 5-1,  and  J  =  N.-  1 ,  and  obtain 

j-1  j  mi 


A“(c.) 


£|'  'bi,ci 


(A. 24) 


where 


bi(cl ’  =  hj  ‘  (1  +  cjyi)c 


'Vi 


(A. 25) 


,  Here  we  will  vary  c.  (sa,),  which  defines  b.(c,  )  for  each  y.,  and  find  the 

12  i  1  'i 

c  (j  =  2,  J^)  that  minimize  2k  (c  ^  ).  Proceeding  in  the  usual  way  by  setting 

da2 

-r—  =  0:  m  =  2.  J  (A. 26) 

dc  m 

m 
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yields  J  -1  equations  of  the  form  below, 

m 


EB  c  .  s  B  .  A  , ;  m  =  2 
mj  j  m, J  +1 
m 

r-  2 


=  2.  J 


(A. 27) 


w  her  e 


in 

B 


m+j  -2ciy, 


(A. 28) 


v*  m  .  ,,  .  -clYi|  - c  1 V i 

B  T  X  1  y  y  b.  -  (1  +c.y . )  e  e 

m,  J  *  1  i  i  l  l  j 


(A. 29) 


The  programs  used  to  solve  this  sot  of  equations  with  the  procedure  described 
are  given  in  Section  A. 5  below. 

A.**  Program  Description  and  Application 

In  this  section  the  FORTRAN  II  programs  vised  to  fit  the  h (y )  data  and  to 
tabulate  the  results  are  given,  .••long  with  results  of  application  to  a  "typical'1 
h(y)  -  Fig.A.l.  There  are  two  programs:  MAGLSQ.  FT  and  MAGTAB.  FT. 

The  first  reads  in  the  h (y )  data  and  obtains  the  least  square  fit.  The  second 

gives  the  constants  a  and  g  ,  and  tabulates  »ho  calculated  results  for  h(y)  and 

J  J 

R(y  >. 

MAGLSQ.  FT  is  listed  in  Fig.  A. 3,  with  input  and  output  in  Tables  A. la  and  A. lb. 


Statements 


10-15 


Description 

Read  in  tabular  array  of  y (j),  h(j)  for  NPT 
points  (Table  A.  la).  This  program  is  re¬ 
turned  to  by  MAGTAB,  hence  it  is  necessary 
to  decide  initially  (K)  whether  to  input  data 
on  start  new  iteration. 

Choose  number  of  max  pt  to  be  fitted  Np(fit  al¬ 
ways  starts  at  j  =  1  >  and  total  number  Nj  of 
a-  coefficients  (top  Table  A.  lb).  Set  up  an  array 
P(j,k)  needed  later. 
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O'  U 


1 


{ 


C  HAGLSG .  FT  J  HT  TO  H(Y)  ^  OR  HAG  ANALYZER.  1  /9/7B. 

C  HAX  VAL  HI  IS  91  VAL  '  S  CIK. OOOOl  NEGLECTE P ( SEE  21  E fC  PLLUW). 

COHHON  NET  .Y.H.NP.NI.JH.C.P.P 

01 HENS I ON  C<9>.H<30>.Y<30>.P<30»16>.P<B,9> 

WRITE <1.1 13)  H 

1  READ< 1 » 1 14)K 

IE(K>3. 10. 13 
3  READ< 1 . 100>NPT 

WRITE(l.lll) 

00  5  J-l .NET 
WRITE ( 1 . 101 ) J. 

5  REAO< 1 . 102)Y< J).H(J) 

10  REAO< 1 . 1 10)NE 

REAO< 1 . 1 12)NI 
JH-NI -1 
H«2*JH 

00  13  J* 1 . NET 
00  12  K-l.H 
P< J.K)-Y< J)**K 
CONTINUE 
WRITE < 1.103) 

REAO< 1 . 106)L 
IE (L  > 1 . 17. 20 
17  REAO< 1 . 107)C0.C0.Cr 

WRITE! 1 . 100) 

N- 0 

19  C1-COTELOAT(N)*CO 
GO  TO  21 

20  REAO< 1 . 103) Cl 

21  z-ioo.»<ci+.ooooov> 

I-Z 

Z-Z-EL0A1 < I > 

H-1000.4Z 

C1-<EL0AT< I )  +  EL  OAT  <  H  > / 1 000 • >/100. 

C<  1  >-Cl 
0-2. *C1 

C  C  <  J  > -A  <  J4 1 >  IN  WRITE-UE.  NOW  EINO  CONSTb. 

00  30  H-2.JH 
PH-O. 

00  23  I-l.NE 
Z-C 1 • Y <  I  ) 

E-0 

IE  <  Z-86 . >22.22.23 

22  E-EXE(-Z) 

23  PH-PH+ <M< I )-< 1 . +Z>»£ >*E*P< 1 »H) 

P  <  H . JH4 1 ) -PH 

00  27  J-2.JH 
PH-O. 

00  23  I-l.NE 
Z-0*Y< 1 ) 

E-0 

IE  <  Z-06 . >24,24.25 

24  E-EXE< -Z ) 

25  PH-PH4E4E! I »H4J> 

27  P  <  H»  J) -PH 

30  CONTINUE 

C  ALL  CONSTS  EOUNO.  SOLVE  EO  ARRAY. 

Fig. A.  i  MAGLSQ.  FT  Lmting 
(cont.  nrxt  page*) 


n  o 


6  l 


PO  70  K-2.JM 

E-B!K»K> 

I*K 

DO  62  N-K.JM 
IF!B<H. J)-E>62. 62.61 
I-M 

E«B!H. J) 

CONTINUE 

t'O  64  JtR.NI 

E-B<K. J) 

p!  R • J)-B! J . J) 

B(I» J)-t 

t'O  66  M-2.JN 

C(M)-R!H.R) 

PO  69  J-K.NI 
B!K,J>-P!K.J>/C!K> 

HO  6C  N-’2»  JH 
ir !M-K )  67  • 68 • 6? 

B<H.J)-B<N.J>-C<H>*R!K*J> 

CONTINUC 

CONTINUE 

CONTINUE 

no  7i  j»2>jn 

C< J)»B< J. JM4  1 ) 

SOLN  CONFUTE  KUN  C<  J>  .  J-2. JM.  WITH  C1-A2  AS  1  NF'UT  »  ANP  Al-1. 

CALC  VAuUf  OK  Nf AN  SR.  PEV. 

P-0. 

PO  75  I  -  1  .  NF- 
Z-C 1  *Y !  I  > 

E-0 

IE(Z-86. 172.72.73 

E»EXF!-Z>  Fill, A. i  (coni.)  MAGLJjQ.  FT  U»t»ng 

BN-0 . 

DO  74  J-2.JN 
BM-BM+C! J)*F( I . J ) 

BM-C*BM  M! 1 >4! l.+Z)*C 
P-B4PH**2 

F  “S0RT<P/FL0AT!  NF’ )  > 

IF<L  > 16.82.00 
WRITE! 1*1 04 ) E 
00  TO  84 

WRITE! 1 . 109>C1  »l  .  P 
N-NK1 

ir<Cl-CF  >19.16.16 
C Ai t  CHA 1 N (  NAG TAB') 

F  ORNAT t ' NO .  KTS  TO  T AB . . NF T - ' I  3 > 

K  ORNA  F  < 13 > 

FORMAT (2X.K5.2.2X.F6. 4) 

FORMAT! ' A2- 'F 10. 6> 

F  ORNAT (  '  H. SO. PE V.- 'FI  0.7) 

F ORNAT ( // ' SI T  L'-IIOM  TO  START  NEW  riT)  ITERATE  A2»' 

'  CHOOSE  A2  ANP  TABULATE.'/) 

FORNAT<  L - ' 1 3 ) 

F  ORNAT  (  '  A2  N1N-  F10.6./'PEL  A2-  F10.6./ '  Ai?  NAX-'F10.6) 
FORNAT<//'NEAN  SQUARE  pro.  IS  FOR  NR  RTS'/'  A2  H.SG.BEV.' 

'  SUN  SO.  PEV. ' > 

FORMA T ( F7. 4. 7X.F I  0. 7. 2X.K 12. V) 

FORMAT! 'NAX  RT  TO  FIT.  NR»'I3> 

FORNAT  < // '  J  Y ( J )  H(J>'/'III  XX. XX  X.XXXX') 

F  ORNAT  <  'TOT •  NO.  A(J)  COEFFICIENTS.  N1-'IJ> 

FORMAT! //'SCT  K--1KM1  TCl  INRllT  FIT  PATAI  NR.NI  VALUES) '/ 

'  OR  TO  START  ITERATION.') 

FORMAT! 'K« ' 1 3) 

END 


> '  F  1 0 . 6  ) 


1  14 
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Ft  HAGLSC 


SET  K--U0J1  TO  INFUT  FIT  DATA!  NP,NI  VALUES1 
OF  TO  STAPT  ITEFATION. 

K  ■  - 1 


0. 

PTS  TO 

TAB. » NPT 

J 

Y(J) 

H<  J> 

1 1 

XX. XX 

x.xxxx 

1 

0.50 

.9600 

2 

1  .00 

.7100 

3 

1 .50 

.2650 

A 

2.00 

.0600 

5 

2.50 

.0055 

6 

3.00 

- .  0 1  30 

7 

3.50 

-.0155 

8 

A. 00 

-.0160 

9 

A. 50 

-.0155 

10 

5.00 

-.01 30 

1  1 

5.50 

-.0100 

12 

6.ee 

-.0080 

1  3 

6.50 

-.0060 

1  A 

7.00 

- • 00A5 

1  S 

7.50 

-.0035 

16 

8.00 

-.0030 

1  7 

8.50 

-.0025 

18 

9.00 

-.0020 

19 

9.50 

-.0018 

20 

10.00 

-.0017 

21 

10.50 

-.0016 

22 

1  1  .00 

-.0015 

23 

11.50 

-.001  A 

2A 

12.00 

-.001 3 

25 

12.50 

-.0012 

26 

1  3.00 

-.0011 

27 

13.50 

-.0010 

28 

1  A. 00 

-.0009 

29 

14.50 

-.0008 

30 

15.00 

-.0007 

T  able  A.  la.  h(y)  Input  Data  for  MAGLSQ.  FT 
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MAX  PT  TO  FIT*  NP-3C 

TOT.  NO.  A<J>  COEFFICIENTS*  NI-9 


SET  L--1J0M  TO  STAFT  NEW  FIT! 


ITEPATE  A2I  CHOOSE  A2  AND  TAFULATE 


L  "0 

A  2  MIN-2. I  see 
DEL  A2-0.020P 
A  2  MAX-2. I 9P0 


MEAN  SQUAPE  DEV •  IS 
A2  N.  SC. DEV. 

2.1500  e. 0009633 
2.17PP  0 • 0009*63 

2.1900  0.0011362 


Note:  A  separate  broad  range  scan  on 
a>  (not  shown)  between  1.0  and 
4.0  found  the  region  2.15  to  2.19 
to  contain  the  minimum.  Shown 
FOR  NP  pt s  here  is  the  scan  of  that  region. 

SUM  SO.  DEV. 

0.000027836 

0.000026866 

0.000038726 


SET  L-1J8JI  TO  START  NEW  FIT/’  ITERATE  A2I  CHOOSE  A 2  AND  TABULATE 


L*  1 

A  2- 2 . 1 700 

M. SO. DEV.-  0.C009A63 


Table  A.  lb.  Kit  Results  lor  MAGl*SQ.  F  I, 


N  =9.  N  30. 

l  p 


57 


16  + 


17  to  19 
20 

21  to  30 
30+  to  71 
71+  to  76+ 

77 


Decide  (L)  whether  it  is  desired  to  iterate  a? 

in  order  to  tabulate  V -^(a^l/N^,  ur  choose 

a  particular  value  of  a^  and  eventually  proceed 
to  MAGTA3.  FT  to  tabulate  results  of  fit. 

Read  in  a,  data  for  iteration  (if  L  0)  (top  Table  A.  lh). 

Read  in  a,  to  obtain  coefficients  and  tabulate  re¬ 
sults  (L  fl  (bottom  Table  A.  lb). 

Find  13  coefficients  (eliminating  roundoff), 
mj 

Solve  array  for  a.  (c.  in  program) 

Calculate  mean  square  deviation  ^(a^)  and  dis¬ 
play  depending  on  L.  value. 

Proceed  to  MAGTAB.  Fl  (only  for  L  1  ) 


MAGTAB.  Fl  is  listed  in  Fig. A. 4, with  output  in  Table  A,  2. 

Stat  ement  Description 

78-  to  89  Calculate  and  tabulate  a  and  g.  (top  Table  A. 2). 

1  J 

89+  to  95  +  Calculate  and  tabulate  h ( y )  and  R(y)  (bottom 

Table  A. 2.  First  throe  columns  are  input  data. 

Table  A. la. Fourth  column  is  the  calculated 
value  of  h(y),  and  fifth  is  h--h(y).  Last  column 
gives  P(y).  Final  operation  in  program  returns 
to  MAGLSQ.  FT  for  further  fitting. 

A. 6  Conclusions 

As  shown  above,  for  ‘he  particular  h(y  >  data  used  an  rms  difference  of 
only  .0009  can  be  obtained  with  nine  coefficients.  This  is  quite  small  com¬ 
pared  to  the  accuracy  with  which  h(y)  can  be  measured.  For  four  coefficients 
the  result  was  .007  3  and  for  seven  it  was  .0019.  Much  was  gained  going  from 
4  to  7,  but  little  thereafter.  The  Nj  9  fit  is  better  than  the  im  asurement  accuracy. 

The  h( y >  fit  is  definitely  sufficiently  accurate  to  justify  its  use  for  the  h(y) 
data  used.  The  g]  coefficients  could  be  employed  directly  in  (A. 18)  (with  K  a  , ) 
to  find  R(y).  Eqs.  (A. 12),  (A. I  3),  (A.  1  5)  and  (A.  16)  then  yield  c  ,  x  and  x(y) 
for  arbitrary  c.  The  only  limitation  on  e.  is  that  it  be  sufficiently  small  to 
justify  the  assumption  'hat  h  depends  only  on  y. 
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r.  r. 


C  MAGT  Afi .  r  T  :  TABULATE  FIT  RESULTS  FROM  MAGLSO.  1/9/78. 

COMMON  NFT.Y.H.NF.NI.JM.C.P 
PIMCNSION  L!V) ,H! 30) . Y130) »f! 30. 16) »G!9> 

C  SET  C!J41>-C!J>.  SO  C( J ) “A <  J ) .  fAK.  COEFT'S  C!J>  AND  G(J). 

DO  78  J-  1  . JM 

78  C!NI+1-J)-C!NI-J> 

C!  1 >-l  . 

TiO  88  JM  .NI 

I-J-2 

Z-l  . 

IF(I)80. 80.79 

79  PO  80  KM  •  1 
Z-Z*<FLOAT<K  )♦  1  .  ) 

80  CONTINUE 
G< J)-Z 

x-o. 

PO  86  JS-J.NI 
I = JS-2 
7-1  . 

IF! I )82.87.ftl 

81  DO  82  K-I . I 

z-zt  t  Float  <  k ) ♦ i . > 

82  CONTINUE 

X-X4Z*C( JS)/C( 2)** < JS-J) 

86  CONTINUE 

G<  J>-X/G<  J) 

88  CONTINUE 
WRITE! 1 . 1 10) 

PO  89  J»1 ,NI 

89  UR  I TE  < 1 • 1 06  >  J. C ! J  > • G I J  > 

coErriciFNTs  tapulatfp. 

TAP.  J. Y ! J) *  M<  J ) ACTUAL »  H  (  J )L ALL .  P1FF  *  ANP  R(Y). 

PM-O. 

WRITE! 1.107) 

PO  95  Jr  1 • NPT 
X»1  . 

Z*G  < 1 > 

PO  90  K-2.NI 
P-F! J.K-t ) 

X-X4P*C!K  > 

90  Z-7+P9G  <  K ) 

P-0 

0*C  (  2  )  *  Y(  J  ) 

IF !0~86.  )93.93.94 

93  P-EXP(-O) 

94  X-Xtp 
7-PA7/C ! 2  > 

P-M< J)-X 
PM-PHt  P**2 

95  WRITE! 1 »I08>J»Y!J) »H!  J) • X. D»  Z 
PM-SORT ( PM/n UAT (NET ) ) 

Z-G  < 1 )/C<2) 

WRITE! 1 . 109) PM. I 
CALL  CHAIN! 'HAGLSQ' > 

106  F0RMAT(I3.2!2X.EI4,6)  ) 

107  FORMAT!//'  J  Y  HACT  MCALC  DIE  F  R! Y ) ' ) 

1 08  FORMAT!  13. 2X.F1>.  2.  J  CM.F  7. 4  >  .2X.F  8.5  > 

109  F  ORMAT  !  / '  RMS  PIFF  ALL  FTS>"'F9.6.  R !  0 )  •  '  F8 . 5  •  //// > 

110  FORMAT!//'  J  A ! J )  G!J)'> 

END 


Fig. A. 4.  MAC"  at  .  FT  Listing 


S9 


J  A<  J  )  G<  J  > 

]  0.100000£^01  0 • 25559  3E^  0 1 

2  0.21  7000F.'»01  0 . 3376  36E*  0  1 

3  -  0 . 602658F.*  0  t  0.13089IE>ei 

4  0 • 353369E*  02  0.530600E+01 

5  -0.456640E*02  -0.1 629  J  8E*  02 

6  0.258635E*02  0.127476E*02 

7  -0. 748371E*01  -0 . 474359E* 0 1 

8  0 . 1 08  7  36E*  0 1  0 • 8494  36E*  00 

9  -0 . 645  356E-0 1  -  0 . 645  356F.- 0  l 


J 

Y 

HACT 

HCALC 

DIFF 

P(Y) 

1 

o.se 

0.9600 

0.9600 

0.0000 

0.70803 

2 

1  .00 

0.7100 

0.7101 

-0.000 1 

0.26541 

3 

1  .50 

0.2650 

0.2648 

0.0002 

0.02889 

4 

2.00 

0.0600 

0.0605 

-0.0005 

-0.04288 

5 

2.50 

0.0055 

0.0044 

0.  CC  1  1 

-0.05594 

6 

3.CC 

-0.0130 

-0.0113 

-0.0017 

-0.05352 

7 

3.50 

-0.0155 

-0.0166 

0.0011 

-0.04628 

8 

4.00 

-0.0160 

-0.0167 

P.0007 

-0.03780 

9 

4.50 

-0.0155 

-0.ei46 

-0.0009 

-0.02992 

10 

5.00 

-0.01 30 

-0.0122 

-0.0000 

-0.02322 

1  1 

s.se 

-0.0100 

-0.0102 

0.0002 

-0.01 762 

1  2 

6.00 

-o.oose 

-0.0084 

0.0004 

-0.01296 

1  3 

6.50 

-0.0060 

-0.0067 

0.0007 

-0.00918 

1  4 

7.0e 

-0.0045 

-0.0051 

0.0006 

-0.00624 

1  5 

7.50 

-0.0035 

-0.0037 

0.0002 

-0.00406 

1  6 

8.00 

-0.0030 

-0.0025 

-0.0005 

-0.00254 

1  7 

8.50 

-0.0025 

-0.0016 

-0.0009 

-0.00153 

18 

9.00 

-0 • 0020 

-0.0010 

-0.0010 

-e. 00089 

1  9 

9.50 

-0.0018 

-0.0006 

-0.0012 

-e. ooosfl 

20 

10.00 

-0.0017 

-0.0003 

-0.0014 

-0.00028 

21 

10.50 

-0.0016 

-0.0002 

-0.0014 

-0.00015 

22 

11.00 

-0.0015 

-0.0001 

-0.0014 

-0.00008 

23 

11.50 

-0.00  14 

-0.0001 

-0.001  3 

-0.00004 

24 

12.00 

-e.coi  3 

-0.0000 

-0.0013 

-0.00002 

25 

12.50 

-0.0012 

-0.0000 

-0.0012 

-e.ceeoi 

26 

1  3.00 

-0.001  1 

-0.0000 

-0.0011 

-e.eoeoo 

27 

1  3.50 

-0.0010 

-0.0000 

-0.0010 

-0.00000 

28 

14.00 

-0.0009 

-0.0000 

-0.0009 

-0.00000 

29 

14.50 

-0.0008 

-0.0000 

-0.0008 

-e.ooejo 

30 

15.00 

-0.0007 

-0.0000 

-0.0007 

-0.000.00 

PKS 

DIFF  ALL 

PT  $■  0. 

000946 

F<  0  )  •  1 

. 1 7785 

Table  A.2.Tabulation  by  MAGTAR.  FT.  N.^9,  N  =30. 

1  P 
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